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Lattice QCD allows us to simulate QCD at non-zero temperature and/or densi- 
ties. Such equilibrium thermodynamics calculations are relevant to the physics 
of relativistic heavy-ion collisions. I give a brief review of the field with empha- 
sis on our work. 



1. Introduction 

High temperature hadronic matter was certainly present in the early uni- 
verse. Relativistic heavy-ion colliders such as RHIC and in future the 
LHC with heavy-ions, produce hot hadronic matter. Lower energy ma- 
chines can produce hot hadronic matter with an appreciable baryon/quark- 
number density - hot nuclear matter. At high enough temperatures this 
hadronic matter is expected to become quark/gluon matter - the quark- 
gluon plasma. Preliminary evidence for the quark-gluon plasma has been 
reported from RHIC and CERN. 

At high baryon/quark- number density and low temperature such as 
might exist in the cores of neutron stars, exotic states of matter, such as 
colour superconducting phases have been proposed. 

In Quantum Chromodynamics (QCD), the accepted theory of hadrons 
and their strong interactions, the most interesting finite-temperature prop- 
erties are non-perturbative. Lattice QCD simulations are the only system- 
atic approach to the study of such phenomena. 

Fig. 1 shows a simplified version of the proposed phase diagram. (For 
an introduction to the phase structure of QCD and predictions from lattice 
QCD, see recent reviews 1 ' 2 ). 

In Lattice QCD, we only really know how to calculate equilibrium ther- 
modynamics. Relativistic heavy-ion collisions are clearly not equilibrium 
processes, but knowledge of equilibrium thermodynamics and the conse- 
quent extraction of thermodynamic quantities and the prediction of the 
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Fig. 1. Simplified phase diagram for QCD at finite temperature and quark- number 
chemical potential fi. 

equation-of-state give us some idea of the type of behaviour we might ex- 
pect. In addition, it can provide input for hydrodynamic models of such 
collisions. 

QCD at finite temperature and/or densities yields information on the 
dynamics of confinement and chiral symmetry breaking and thus enhances 
our understanding of QCD. 

2. Lattice QCD at finite temperature 

QCD is quantized by functional integral methods. To make these integrals 
well-defined, we rotate to imaginary time - Euclidean space-time. The quark 
fields are then defined on the sites of a 4-dimensional hypercubic lattice, 
the gauge fields = exp^gA^) on its links. 



(1) 
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where the action S = Sf + S g . The simplest lattice implementation (which 
preserves gauge invariance) has 

S s= 0^(l-^ReTr n UUUU), (2) 

where j3 = 6/g 2 , and Sf = Y^ites m )V>- 

Simulations are performed by replacing the fermion fields with boson 
fields (pseudo-fermions) which results in replacing the Dirac operator by its 
inverse. Adding a 'kinetic' term for the gauge fields, allowing them to evolve 
in a fictitious 'time', turns Z into the partition function for a system of clas- 
sical particles. One then uses molecular-dynamics techniques to effectively 
'evaluate' the integrals, by numerically integrating the equations of motion 
for this system. Bringing the system in contact with a heat-bath at regular 
intervals assures ergodicity and compensates for the lack of dynamics for 
the pseudo- fermion fields. 

To compensate for the doubling problems for the fermion fields fre- 
quently requires taking fractional powers of the fermion determinant. This 
is performed by using rational approximations to fractional powers of the 
Dirac operator. Use of a global Metropolis accept/reject step removes the 
discretization errors in the numerical integration of the equations of motion. 
(For a recent description of simulation methods with dynamical fermions 
see 3 ). 

If we use a lattice of temporal extent 1/T and a spatial extent much 
larger than this, and demand periodic boundary conditions for the gauge 
fields and antiperiodic boundary conditions for the fermion fields in the time 
direction, Z = Z(T) the partition function for lattice QCD at temperature 
T. 

3. The finite temperature transition 

For massless quarks, chiral symmetry is restored at the finite tempera- 
ture transition. (V>V) below the transition and (ipip) = above the 
transition. Hence this transition is a phase transition. Arguments based on 
dimensional reduction predict that this phase transition is a second order 
phase transition (critical point) for Nf = 2, in the universality class of the 
3-dimcnsional 0(4) spin model. 4 For non-zero quark mass it is expected to 
weaken to a crossover with no real phase transition. For Nf > 2 this tran- 
sition is predicted to be a first-order phase transition which is expected to 
remain first-order for small quark masses, becoming a crossover for larger 
quark masses. 
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In the real world where Nj = 2 + 1 (u,d,s), the important question 
is whether the strange quark mass is light enough for the transition to be 
first-order or whether it is a crossover. Recent lattice simulations indicate 
that the strange quark mass is too large, and the transition is a crossover 5 
(see Fig 2). This extends earlier work of 8 (NB the Nf = 3 point agrees with 
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Fig. 2. The chiral critical line at zero baryon-number density (from Ref. J ). 

our own simulations.) 

Simulations with larger lattices, finer lattice spacings and improved ac- 
tions have recently been reported by. 6 Here the u, d and s quark masses 
are chosen by fixing mx/m^ (and fx/ fn) at their physical values. Making 
several choices of masses restricted in this manner, they were able to extrap- 
olate to the physical quark masses and below. They confirm the observation 
that the physical transition is a crossover and not a phase transition, in sim- 
ulations which probe much closer to the continuum limit than in. 5 The same 
simulations give new estimates for the temperature(s) of this crossover. 7 

Since the strange quark appears too massive to control the nature of 
the transition, it is useful to study the 2-flavour case, to see if the nature 
of the m = phase transition agrees with the above prediction. 

Determining the nature of the 2-flavour phase transition has proved 
difficult, since is expected only at m = 0, and finite volume effects appear 
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to be large. 9-12 Standard (including highly improved) lattice actions forbid 
simulations at to = 0, and small mass simulations are very expensive. No 
such simulations have been performed at masses small enough to uncover 
the to = results. 

We have performed Nf = 2 simulations using the xQCD action which 
allows simulations at zero quark mass. 13 Since it is a staggered quark action, 
the reduced flavour symmetry leads us to predict that the phase transition 
should be second-order and in the universality class of the 3-dimensional 
0(2) spin model. To accommodate the finite volume effects, we compare 
our chiral condensate measurements to the magnetization of the 0(2) spin 
model also on finite volumes, as shown in Fig. 3. The agreement is excellent. 
Good agreement is also obtained with correlation lengths and susceptibili- 
ties. 
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Fig. 3. Chiral condensate (points) near the chiral phase transition for lattice QCD, 
fitted to the magnetization (solid line) for the 0(2) spin model (from Ref. 13 ). 
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Fig. 4 shows recent estimates of the transition temperature T c from 
lattice QCD simulations, compared with an experimental estimate from 
RHIC. All lattice simulations are consistent with the chiral and deconfine- 
ment transitions being coincident. 
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Fig. 4. Recent lattice estimates of T c compared with experiment (from Ref. 14 ). 



4. The equation-of-state for hot QCD. 

The equation-of-state (EOS) expresses the pressure p, the entropy density s 
and the energy density e as functions of the temperature T. The free energy 
density, pressure and entropy are given by 

/ = ~lnZ(T), p=-f ,s=^, (3) 
respectively. The energy density is not an independent quantity but is given 

by 

e = Ts- p. (4) 
Z(T) is calculated on the lattice by numerically integrating 

dlnZ I Q \ M 
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To obtain T as a function of requires knowing the lattice spacing as a 
function of (3. This can be obtained by measuring physical quantities as 
functions of (3 at zero temperature. 

Knowledge of the EOS is needed as input to models for the evolution 
of the hot hadronic matter in relativistic heavy-ion collisions. Fig. 5 is a 
graph of these quantities from the MILC collaboration. 15 
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Fig. 5. The EOS from Lattice QCD simulations (MILC preliminary: status Lat- 
tice2006 16 ). I = e - 3p 



For earlier work on the finite temperature EOS see for example. 17 Karsch 
et al. have shown that the low temperature behaviour is well modeled as 
a non- interacting gas of hadron resonances. 18 See also 19 for other work on 
the QCD EOS. 

5. Meson spectral functions at finite T 

Meson spectral functions yield information about the propagation of 
mesons, or excitations with mesonic quantum numbers in hot hadronic 
matter 20 (This review gives references to earlier works). Even just above 
T Cl the quark-gluon plasma is a strongly interacting fluid and mesonic states 
survive. Not only do the spectral functions have information about hadronic 
stability at high temperatures, but they also have information about trans- 
port coefficients and dilepton production. At zero temperature, the spectral 
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function is just the momentum space propagator. The Euclidean-time me- 
son propagator G(r,p,T) is related to the spectral function a(uj 7 p,T) by 



The main difficulty has been that, since the spatial extent of the lat- 
tice must be much greater than its temporal extent, the temporal extent of 
the lattice is typically ~ 10 or less in lattice units. This gives a rather poor 
estimate of a. Recently simulations have been performed on anisotropic lat- 
tices on which the spatial lattice spacing is much greater than the temporal 
lattice spacing. 21 ' 22 This allows for many more points in the time direction, 
while still keeping the spatial extent of the lattice much greater than its 
temporal extent. This leads to better estimates for the spectral functions. 

The charmonium spectral functions have been of particular interest. It 
had been suggested that one signal for the quark-gluon plasma phase could 
be the 'melting' of charmonium. This would reveal itself as a drop in the 
charmonium production relative to the production of DD pairs and similar 
states. 

Current simulations by the TrinLat collaboration using anisotropic lat- 
tices with dynamical light quarks, indicate that the J/ip and rj c survive 
above T c and melt somewhere between 1.3T C and 2T C . The %c states ap- 
pear to melt at or below 1.3T C (see Fig. 6). 

At the LHC, even higher temperatures are expected, and it has been sug- 
gested that these might be high enough to melt bottomonia. 23 We plan to 
measure NRQCD bottomonium propagators (and hence spectral functions) 
on the TrinLat configurations to determine these melting temperatures. 

6. Transport coefficients 

To use hydrodynamic models for hadronic matter in relativistic heavy-ion 
collisions requires knowledge of the transport coefficients, shear viscosity 
r], volume viscosity £ and thermal conductivity k. These viscosities are 
expressed in terms of Green's functions of the stress-energy tensor. 24 ' 25 



These real-time Green's functions are obtained from their lattice (Eu- 
clidean time) counterparts through the spectral function (j(uj). Quenched 
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Fig. 6. Spectral function p = cr for the J/ip at T c , 1.3T C and 2T C (from Ref. 22 ). 

lattice results have been obtained for these quantities. 26 In terms of the 
spectral function, r\ — tt lim^^o 

This means that one needs the spectral function near lo — 0, where it 
is least well known. Thus the determination is difficult, and systematic as 
well as statistical errors are large as shown in Fig. 7. The same paper finds 
(k0. Clearly there is a long way to go. 

More recently, improved methods have been used to calculate the elec- 
trical conductivity of the quark-gluon plasma. 27 

7. Lattice QCD at finite baryon number density 

Finite baryon/quark- number density is implemented by introduction of a 
quark-number chemical potential fi. On the lattice this is achieved by mul- 
tiplying each of the links in the +t direction by e M and each of those in the 
—t direction by e _Al in the quark action Sf. 

Integrating out the fermion fields gives the determinant of the Dirac 
operator which is complex, with a real part of indefinite sign. Standard 
simulation methods, which rely on importance sampling, fail for such sys- 
tems. 

Some progress has been made in circumventing these problems for small 
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Fig. 7. r//s as a function of temperature. Dashed line is at 1/4.7T (from Ref. 26 ). 



/is close to the finite temperature transition. Methods for doing this fall into 
several classes. 



• Analytic continuation: In the simplest case people simulate at imag- 
inary jits, where the fermion determinant is real and positive. The 
results are fitted to a power series in /i 2 , which allows continua- 
tion to real /i. 28 ' 29 Fancier analytic continuation methods have also 
been used. 30 

• Power series expansions: These are similar in spirit to the analytic 
continuation methods. The exponential of the action is expanded in 
powers of /i 2 . The coefficients of this expansion are then observables 
whose expectation values are measured in \i = simulations. 31,32 

• Reweighting methods: These start from a quark action with a pos- 
itive fermion determinant, and reweight measurements by the ratio 
of the original fermion determinant to this positive fermion deter- 
minant. One then divides by the expectation value of this ratio of 
determinants. 33 

• Phase quenched methods: One simulates using the magnitude of 
the fermion determinant ignoring the phase. For small enough \i 
on a finite lattice, the phase is small enough that this should yield 
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the same phase structure as the full simulation. 3 
• Canonical ensemble methods: The fermion determinant is projected 
on to states of fixed quark number. One deals with any sign prob- 
lems by reweighting. 35 

For 3-flavour QCD, it has been argued that the critical point observed 
at zero chemical potentials, where the transition weakens from a first-order 
phase transition to a crossover, would move to larger quark masses as ji is 
increased from zero, becoming the sort-after critical endpoint. Recent work 
using analytic continuation from imaginary /i (de Forcrand and Philipsen 28 ) 
and simulations in the phase-quenched theory (Kogut and Sinclair), 34 in- 
dicate that this does not happen. Instead, the critical mass appears to 
decrease with increasing /x. 

Fodor and Katz, using the reweighting method, claim to find the criti- 
cal endpoint, for physical quark masses. 33 Their estimate of its position is 
T = 162(2) MeV and fi = 120(13) MeV. Since (J, > m n /2, this is beyond the 
reach of analytic continuation, phase quenched and series expansion meth- 
ods. It remains to develop other methods which could check this. Fig. 8 
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Fig. 8. Binder cumulant for the chiral condensate as a function of /ij (from Ref. 34 ). 
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shows the behaviour of the Binder cumulant in the phase-quenched theory. 
If there were a critical endpoint at (small) finite fxj , this Binder cumulant 
would decrease through the Ising value 1.604(1) (dashed line). Instead it 
increases indicating that there is no critical endpoint in this range of 



8. The EOS at non-zero T and /x 

At finite temperature and small /x, the series expansions in terms of \x also 
enable one to calculate such quantities as p, s and e and hence to study the 
equation-of-state in terms of T and /z. 36 Similar results have been obtained 
using multi-parameter reweighting techniques. 37 

The pressure p can be obtained from measurement of the quark-number 
density p, since 

< 8 > 

The pressure at zero p is calculated as above. That at finite fi can be 
obtained by integrating the previous equation, e can also be calculated, but 
requires knowledge of the running of the coupling constant and quark mass. 
Results from the Bielefeld group are shown in Fig. 9 



9. Summary and Conclusions 

Lattice QCD at finite temperature can probe the nature of the phase tran- 
sition from hadronic matter to a quark-gluon plasma. For physical quark 
masses this appears to be a crossover without a true phase transition, in- 
fluenced by the second-order transition at m u = = 0. 

The QCD equation-of-state and other equilibrium thermodynamic 
quantities measured in lattice QCD simulations provide input for an under- 
standing of the non-equilibrium thermodynamics at RHIC and the LHC. 
Finite temperature (lattice) QCD probes QCD dynamics such as confine- 
ment and chiral-symmetry breaking. 

Charmonium spectral functions show that the J ftp remains intact at 
the finite temperature transition T c , but dissociates below 2T C . This should 
produce a reduction in J ftp's from the most energetic processes at RHIC. 
Bottomonium spectral functions should show similar behaviour at the even 
higher temperatures of the LHC heavy-ion program. 

Just above T c , the 'deconfined' phase is a strongly-interacting fluid with 
low viscosity. Transport coefficients (including viscosity) need to be calcu- 
lated as input for hydrodynamic models. Early attempts have been made 
in lattice QCD simulations, but such measurements are difficult. 
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Sign problems hamper simulations at finite baryon/quark- number den- 
sity. Some progress has been made for small /is close to the finite temper- 
ature transition. The equation-of-state has been determined in this high- 
temperature low baryon-number-density regime. More work is needed to 
observe the critical endpoint, which is the most striking feature expected 
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in this region of the QCD phase diagram. This is a region accessible exper- 
imentally to lower-energy relativistic heavy-ion collisions. 

New methods will be needed if one is ever to reach the high baryon- 
number densities needed to understand the physics of neutron stars. At such 
densities one expects such exotic states of matter as colour superconductors. 
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